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0 Method and apparatus for storing, accessing, and processing voxehbased data. 



® The present invention concerns a method of and apparatus for storing to and accessing from a memory 
device, one or more voxels of a beam disposed along one or more of a plurality of storage/retrieval directions of 
a 3-D matrix array and 3-D discrete voxel space. Voxels are stored into memory space by mapping the voxels 
atong one of the plurality of storage/retrieval directions, into a plurality of independently successful ntemory 
modules in memc^y storage space. Each memory module is indexed by a memory module index and has 
internal memory cell addresses. The mapping is carried out in accordance with a linear skewing function which 
expressed in terms of x, y, and z coordinate directions and integer n. The voxels can be retrieved along one or 
more of a plurality of storage/retrieval directions by de-mapping one or more of the voxels from the memory 
storage space into 3-0 voxel space using spadal parameters and integer n. The demapping operations of the 
present method are also used in order to determine the voxel depth measures of voxels along one of the 
■storage/retrieval directions. ^ > ^ i 
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METHOD AND APPARATUS FOR STORINa ACCESSING. AND PROCESSING VOXEL-BASED DATA 



Reld of the Invention 

The present invention relates generally to methods of and apparatus for storing in and accessing from a 
memory device, three-dimensional data. More particularly, the present invention relates to methods of and 
5 apparatus for storing and accessing voxel-based data of voxel-based imagery from a three-dimensional 
frame buffer memory. 



Setting for the Invention 

10 

The need to model objects and processes in two-dimensional and three-dimensional (hereinafter 2-D 
and 3-D, respectively) space is manifest, and the advantages of doing so using computer graphic systems 
is equally as clear. 

Presently, a variety of computer graphic systems are known, and for the greater part, have been 
IS characterized in terms of the data they manipulate. For example, there are systems that manipulate data 

represented on the continuous object (geometric) space R^. there are those that operate on the discrete 

pb(el-image plane Z*. and there are those that work on the discrete voxel-image space 2?. 

Thus, as a result of advances in computational science and technology, it has recently become quite 

conventional to use three-dimensional computer graphic systems having surface and sofid modelling 
20 capabilities. As with other types of graphics systems, such types of systems can be characterized by the 

fact that the 3-D object or scene is represented by a geometric model in 3-D continuous geometric space. 
Such 3-D geometric models are typfcally constructed from continuous 3-D geometric representattons. 

Including for example, 3-D straight fine segments, planar polygons, polyhedra. cubic polynomial curves. 

surfaces, and volumes, circles, and quadratic objects such as spheres, cones, and cylinders. These 3-D 
25 geometrical representations are used to model various parts of an object, and are expressed in the form of 

mathematical functions evaluated over particular values of coordinates in continuous space. Typically. 

the 3-0 geometrical representalions of the 3-D geometric model are stored in the fomiat of a graphical 

display list. 

Once modelled using continuous 3-D geometrical representations, the object can be graphically 

30 displayed along a particular view direction by repeatedly scan-converting the graphical display list 

However, such 3-D geometrical-based graphic systems are not without shortcomings and drawbacks. 
For example, when viewing the 3-D geometrical model of the 3-0 object, only the surfaces visible to an 
observer ateng a particular viewing direction shouU be displayed. However, to ensure this, the hidden 
surfaces must be detected using surface detection algorithms, and removed by a computationally intensive 

35 process refen-ed to as "hidden surface removal". Also, for every small modification, such as for example, a 
change in viewing direction, the graphical display list must be manipulated and repeatedly scan-converted, 
which is both time consuming and computationally intensive, further limiting the technk^ue. 

Another major drawback.of 3-D geometrical-based graphics systems is that projection of a 2-D view of 
the 3-D geometrical model, and rendering and manipulation thereof, are dependent on scene complexity. 

40 In addition, in some applications, like in medical imaging, biology, geology and 3-D image processing, 
the data is already in a voxel-based form and no geometric information is available. 

As an alternative to 3-D geometrical-based systems, several 3-D voxel-image space systems have been 
developed. In general, these voxeWmage space, or voxel-based, systems are based on the concept that a 
3-0 inherently continuous object or scene is discretized. sampled, and stored in a large 3-D Cubic Frame 

45 Buffer (CFB) memory comprising unit cubic ceUs called volume elements, or "voxels*. Three-dimensional 
obiects are digitized and a regulariy spaced 3-D matrix array of values is obtained and stored in the Cubic 
Frame Buffer. The 3-D array of voxels in the Cubic Frame Buffer is loaded with experimental date obtained 
by a 3-D scanner such as a CAT or MRI medical scanner, or from a voxel date base. Alternatively, the 
digitization is perfomied by scan-convertlng a 3-D geometrical modal. Both experimental and synthetic 

50 samplings yield a cellular tessellation of the original volumes, and the resulting cellular cubic memory 
provides a real yet discrete model of the 3-D scene or object. 

Consequently, the voxel representetion, which is fundamental to all 3-D voxel-based systems, is very 
effective for the traditional ^pBcations employing empirical 3-D voxel Imagery, such as medicine (e.g. 
computed tomography and ultrasounding), geology, btology, and 3-D image processing {e.g. time varying 2- 
0 images). 

2 
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The advantages of using voxel representation of objects and scenes are numerous, even if used solely 
for synthetic graphics, i.e. 3-0 scan-conversion generated voxel ot)jects. The projection of voxehbased 
images aiong a view direction displays only the visible surfaces and implicitly ignores the hidden surfaces 
of 3-0 voxel-images. Aiso« for every small modification, there is no need of repeatedly scan-converting and 

s manipulating **graphical display lists** as required in 3-D geometrical-based graphics systems. 

Rather, a change can be made directly to the Cut>ic Frame Buffer images. i.e. objects can be directly 
edited by deleting, masking out, or adding voxels with minimal computational effort. Also, when using a 
voxel-based approach, projection, rendering and manipulation are independent of the scene complexity, and 
it is normally unnecessary to employ conventional display lists, repeated scan-conversions, and the 

10 computationally intensive hidden surface removal or surface detection algoritiims. 

To date, there have been several proposals for 3-D voxel-based graphics systems, all of which have 
one thing in common, namely, the utilization of high performance multi-processing architectures in order to 
process tiie huge amounts of voxel-represented data in the cellular cubic memory. Examples of such prior 
art systems, are disclosed in the papers "An Architecture for the Real-Time Display and Manipulation of 

IS Three-Dimensional Objects" by S.M. Goidwasser et al. published in Proc. Int'l. Conference on Parallel 
Processing at pages 269-274. and in "A Generalized Object Display Processor Architecture" published by 
S.M. Goidwasser In I.S.E.E. Computer Graphics and Applications . Volume 4, Number 10, at pages 43-55: in 
the publication •'The Graphics PARCUM (Processing Architecture Based on Cubic Memory) System: A 3-D 
Memory Based Computer Architecture for Processing and Display of Soiid Models" by 0. Jacket, published 

20 in Computer Graphics Forum, Volume 4, Number 4, 1985 at pages 21-32: and in tiie publication entitled '*A 
Three-Dimensional Display Method for Voxel-Based Representation" by T. Ohashi et al. pMJblished in Proc. 
Eurographics '85 . September 1985 on pages 221-232. 

In response to tiie numerous shortcomings and drawbacks of prior art voxel-based systems, applicants 
have proposed an alternative computer graphic system based upon voxel representation of 3-0 objects In a 

25 large 3-D memory which has been described, for example, in tiie following publications: 

•A 3-D Cellular Frame Buffer, *• Arte Kaufman and R. Bakalash. in Proc. EUROGRAPHICS '85 . Nice, France. 
September 1985, pp. 215-220: 

"Memory Organization for a Cubic Frame Buffer," Arie Kaufman, in Proc. EUROGRAPHICS '86 . Lisbon. 
Portugal, August 1^. pp. 93-100; 
30 "Towards a 3-D Graphics Workstation." Arie Kaufman, in Advances in Graphics Hardware t, W. Strasser 
(Ed.), Springer-Verlag. 1987. pp. 17-26; " 

"Voxel-Based Architectures for Three-Dimensional Graphics." Arte Kaufman, in Proc. IRP |^ 10th World 
Computer Congress . Dublin. Ireland, Sept 1986. pp. 361-366; and 

"CUBE - An Archltec^re Based on a 3-D Voxel Map," Arie Kaufman and R. Bakalash. in Theoretical 
35 Foundations of Computer Graphics and CAD . R.A. Earnshaw (Ed.). Springer-Verlag. 1988. pp. 689-701 . 

As dlsdtSed In tiie above publications and generally illustrated in Rgs. 1 and 2, the proposed 3-D 

computer graphic workstation is based upon the 3-0 voxel-based representation of objects within a large 3- 

0 memory referred to hereinafter as a 3-D Cubic Frame Buffer, which comprises specially organized 

memory modules containing a cellular array of unit cubic cells called voxels. The workstation i is a multi- 
40 processor system witii three processors accessing the Cubic Frame Buffer 2 to input sampled or synthetic 

data, to manipulate, to project, and to render tiie 3-D voxel images. 

In general, the processors Include a 3-D Frame Buffer Processor 3, a 3-D Geometry Processor 4, and a 

3-0 Viewing Processor 5. The 3-D Frame Buffer Processor 3 acts as a channel for 3-D voxel-based images 

which have been "scanned** using a 3-D scanner 6 such as CAT and MR! medical scanners. The 3-0 
49 scanned voxel-based images are the primary source of Cubic Frame Buffer data. Once tiie voxel images 

are stored in tiie Cubic Frame Buffer 2. they can be manipulated and transformed by their 3-D Frame 

Buffer Processor 3, which also acts as a monitor for 3-D interaction. 

The 3-D Geometry Processor 4 samples and thereafter converts or maps 3-0 continuous geometrical 

representations of a 3-D object, into their 3-D discrete voxel representation within a Cubic Frame Buffer 2. 
so Notably, the 3-D continuous geometrical representations comprise a set of matitematical functions which, as 

a whole, serve as a 3-D model of the 3-0 object Together, tiie sampling and conversion (i.e. mapping) 

process is typically referred to as a "scart-conversion'* process. 

Referring to Fig. 3, in particular, a general overview of the 2-0 and 3-0 scan<onversion processes is 

given in terms of (I) mapping from continuous 3-D geometi'lc models to 2-D discrete pixel image space, and 
56 (U) mapping from continuous 3-0 geometric models to 3-0 discrete voxel-image space, respectively. In the 

above-described 3-0 voxel-based systems, the 2-0 scan-conversion process illustrated in Rg. 3 is not 

carried out as such prior art processes are strictiy limited to 2-0 image data based generation and 2-0 pixel 

image modelling, whereas in contrast, tiie 3-0 scan-conversion process provides robust 3-D image data 

3 
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base generation and 3-D voxel imago modelling. 

The 3-0 Viewing Processor 5 examines the voxels in the Cubic Frame Butter 2 in the specified view 
direction which can be one of a variety of directions. By taking into consideration depth, translucency and 
color values, the 3-D Viewing Processor 5 generates a 2-D projection (i.e. video pixel image) of the cubic 
5 frame voxel-based image, inputs the same into a conventional 2-D Frame Buffer 7. which in turn is scanned 
by a conventional video processor 8. thereby updating a video screen 9 with the 2-D pixel image. The 
Viev»ing Processor 5 also performs the viewing transformations, hidden surface removal, and allows for 
different viewing angles, or directions. „ „ - ^ ^ .h« 

Producing in real-time 2-D projections of images stored in the Cubic Frame Buffer 2, depends on the 
10 specified view direction, the view position of the observer, and the viewing parameters, and involves 
handling a huge amount of data. In fact, in order to generate an orthographic projection of a moderate 
resolution (i.e. 5t2 x 512 x 8) the data flux through the Viewing Processor 5 is enonnous. 

In the above-described publications, applicants proposed a method whereby the Viewing Processor 5 
would employ a sequence of (Voxel Multiple-Wfrite) processing units which would simultaneously process a 
IS full beam of voxels along a specified viewing direction, (each voxel in a separate processing unit), and 
select the first opaque voxel closest to the assumed observer viewing along that viewing direction, in a time 
which is proportional to the log of the length of the beam. A voxel value defined as a transparent color 
vrauld be used to disable all the processing units holding transparent voxels, all the units holding opaque 
units would be allowed to access a common bus, each attempting to place its opaque voxel value, one byte 
at a time, onto the bus. The opaque value closest to the observer wins and its voxel value is transferred to 

the 2-D Frame Buffer 7. „ ^ « • i 

Thus by retrieving beams of voxels alter beams of voxels from the Cubic Frame Buffer 2 using voxel 
multiple wife bus. a 2-D projection of the Cubic Frame BuHer image could be generated onto a 2-0 graphic 

^^''to "alleviate the memory accessing (i.e. low storagefretrieval speed) problem, applicants proposed in the 
above-referenced publications, a specific 3-D "symmetrically-skewed" memory storage orgwiiation for the 
Cubic Frame Buffer 2. which coukl provide conflict-free access to the Cubic Frame Buffer 2 atong 
storage^etrieval directions parallel to the X. y and z principal axes thereof. 

As illustrated in Rg. 4A. the skewed memory storage organization invoh/ed sectioning the 3-D voxel 
space of the Cubic Frame Buffer into diagonal parallel sub-spaces having a 45 angle with respect to the 
main principal axes planes. According to this proposed memory organization scheme, a voxel having space 
coordinates x y and z is mapped into the k-th memory module of the Cubic Frame Buffer according to the 
skewing function, k = (x + y + z) mod n. where n is the length of one side of the nxnxn cubic memory 
buffer, and where each internal memory cell indtoe (i. j) within each k-th memory module is given by i = x. 
and j = y Since two coordinates are always constant atong any beam of voxels in the Cubic Frame Buffer 
regardless of the viewing direction, the third coordinate guarantees that only one voxel from the beam 
resides in any one of the memory modules. 

According to this proposal, the physical memory of the Cubic Frame Buffer is divided Into n memory 
modules, as illustrated in Rgs. 4B through 4E. To ensure that no two voxels of a beam reside in the same 
40 module, and therefore, that all the voxels of a beam can be retrieved simultaneously 0* in parallel) in one 
memory cycle, each memory module has its own independent address, k. 

Thus using this proposed memory organization, a full beam of n voxels stored away -scrambled in n 
different memory modules. couW be. in theory, simultaneously accessed for all three pnnapal viewing 
directions. i.e. along the £x. ly, tz principal axes. However, several significant problems arise from the 
45 user/viewer's point of view, when using such a memory storage scheme. Specifically, using the proposed 
memory storage scheme, the tme voxel-based image in discrete voxel space Z>, is "scrambled" and stored 
away in n different memory modules each having nxn memory cells, with address indices fi. D thrt differ 
from the physfcal x, y and z coordinates of the voxels. Thus, as a result of the mapping operation of voxe 
beams into computer-storage space N». the user/viewer is wholly incapable of specifying (in terms o 
so physically interpretable measures) that a partfcular beam of voxels along a partteular storage/retneval 
direction, is to be retrieved from the skewed Cubic Frame Buffer. 

In connection with the above-described beam retrieval problem, it is noted that in the above pnor art 
publications, this particular problem was not only not addressed, but also there was no disclMure. teaching 
or suggestion as to how a beam of voxels mapped into and stored within in memory modules of Cutac 
59 Frame Buffer (according to the above^lescribed "skewed" memory storage scheme), might be retrieved 

^^'^Tlwrt^-besS system proposed in applicants' prior art publications, is not without further limitations 
and shortcomings. In particular, the proposed system is Hmited to parallel projection of 3-0 voxel images. 
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along one or more of the six principal view directions: +x. -x, +y, -y. +z. -z orthographic views paraltef to 
the principal axes. Thus, viewing of the 3*0 voxel image in the Cubic France Buffer from more than the six 
principal viewing directions, requires that the 3-0 voxel image in the CulDic Frame Buffer is first transformed 
(i.e. rotated) using the 3-D Frame Buffer Processor (an operation which distorts the 3-D voxel image), and 

5 thereafter viewing it along one of the principal orthographic directions. It is only through such a process that 
arbitrary parallel projections are achievable. 

While the prior art publications have described the Voxel Multiple Write Bus as employing a sequence 
ol processing units which simultaneously process a full beam of voxels and seiect the first opaque one in a 
time which is proportional to the log of the length of the beam, the details for canrying out such a concept 

fo have not been disclosed or known. 

In addition, while prior art publications have sketched a general proposal of how tte Cubic Frame Buffer 
might be organized in order to map voxel values into the memory cells of a plurality of memory modules, 
there has been no disclosure or suggestion as to how the Cubic Frame Buffer might be accessed in order 
to retrieve in a parallel fa^ion, a beam or an interval of voxels from the Cubic Frame Buffer, along a 

15 specified storage/retrieval direction. 

There also has been no disclosure whatsoever of the details of an actual functional embodiment of the 
Cubic Frame Buffer memory storage organization and its addressing system, nor has such a system been 
known. 

Moreover, while the general manner of possible operation of a Cubic Frame Buffer, a Voxel Multiple 
20 Write Bus. and a memory addressing system of a possible voxel-based system has been briefly discussed 
in the prior art publications, there has not been a dieclosure. teaching or suggestion as to how the 
components of an operational Cubic Frame Buffer, Voxel Multiple Write Bus and menwry addressing 
system could be configured to provide an operational 3-D voxel system. 

In view, therefore, of prior art 3-D voxel-based graphic systems, there is a great need for 3-D voxel- 
25 based graphic systems aruj methods which are capable of storing in and retrieving from 3-D skewed 
computer storage space, beams of voxels along a wide variety of storing/retrieving (i.e. viewing) directions 
in 3-D discrete voxel space Z^. requiring specification of the beam in voxel space, and solely In terms of 
spatial parameters having a physical interpretation to the user. 

In addition to 3-0 voxel-based graphic systems which are fast computationally efficient, and organized 
30 to simplify hardware and software implementation, there is also a great need for 3-D voxel-based graphic 
systems and methods capable of generating parallel projections of 3-D voxel images from a various number 
of viewing directions, without distorting or othenwise degrading the quality of infonnation of such 3-D voxel 
images, by requiring transformation (i.e. rotation) of the 3-D voxel image stored in the Cubic Frame Buffer, 
prior to projection. 

35 AcconJingly. it is a principal object of the present Invention to provide a method of storing in and 
retrieving from 3-D computer storage space, beams of voxels parallel to one of several possible conflict-free 
storage/retrieval directions in 3-D discrete voxel space Z^, requiring specification of the beam of voxels in 
Z^ voxel space, solely in terms of spatial parameters having a physical interpretation to the user/observer. 

40 

SUMMARY OF THE INVENTION 



According to the present inventfon. a method is provided for storing and retrieving one or more voxels 
45 of a beam disposed along one or more of a plurality of storage/retrieval directions of a 3-D matrix array and 
3-D discrete voxel space. A 3-D disaete voxel space is specified in terms of x. y and z coordinate 
directions. In general, the method comprises storing the voxels into memory storage space, by mapping the 
voxels along one of the plurality of storage/retrieval directions, into a plurality of independently memory 
modules In the memory storage space. Each memory module is indexed by a memory module index and 
50 has internal memory cell addresses, in a prefenred embodiment, the mapping is carried out in accordance 
with a linear skewing function expressed in terms of the x. y and z coordinate directions and integer n. The 
voxels are retrieved along one of the plurality of storage/retrievable directions in 3-D discrete voxel space, 
by demapping one or more voxels from the memory storage space into the 3-D voxel space using spatial 
parameters and integer n. 

95 in the preferred embodiment the mapping of the voxels along the beam parallel to the storage/retrieval 
direction, is into a plurality of n independently accessible memory modules whfch are indexed with memory 
module indices kaO, 1. ...n-1. Each memory module contains n^ memory cells and each memory cell is 
specified by an internal memory nfK)dule indices. The mapping is carried out in acc(»'dance with the linear 

9 



6NSDOCID: <EP 036a42SA2 I > 



EP 0 368 425 A2 



10 



IS 



skewing furction of the general form k = (ax + by + cz) mod n. where x, y and z of are to coordinate values of 
each voxel in 3-D discrete voxel space, where i='x and j = y. and where a. b and c are integer coefficients 
and n Is prime relative to the integer coefficients a, b and c and linear combinations of the integer 
coefficients a. b and c. Notably, the number of storage/retrieval directions in the n^ matrix array in 3-D 
discrete voxel space is based upon the value of integers a. b. c and n. 

To retrieve the voxels along one of the plurality of storage/retrieval directions in 3-D voxel space, the 
following procedure is performed. For each voxel along the beam parallel to the storage/retrieval directions, 
a voxel depth measure representative of the depth of each voxel In the n^ matrix array is determined on the 
basis of one or more of the x. y and 2 coordinate values of the voxel, the integer coefficients a. b and c and 
the module index k of the memory module Into which the voxel has been mapped. Then, for each voxel 
along the beam parallel to the storage/retrieval direction, the memory cell indices i and j are determined on 
the basis of the voxel depth measure and the x. y and z coordinate values corresponding to each voxel. 
Then, the voxels along the specified storage/retrievable directions are retrieved using the memory cell 
indices i and j determined above and the module index k. 

In one embodiment of the present invention, six storage/retrievable directions corresponding to the 
principal axis +-x, +-y and +-z of the matrix array, are provided by setting the integer coefficients a, b 
and c « 1 in order to deternnlne the voxel depth measures for voxels along the beams parallel to the 
various principal axes, various demapping functions are utilized. Notably these demapping functions are 
derivable from the linear skewing function used in the memory storage mapping operations method. 
20 In the preferred embodiment. 26 storage/retrieval directions of the 3-D matrix array in 3-D voxel space 
are provided, by setting integer coefficiente a = 5. integer coefficient b = 2 and integer coefficient c-1. In 
this case, the voxel depth measures for voxels along beams parallel to any one of the twenty six 
storage/retrieval directions, require particular demapping functions which are derivable using modular 
arithmetic, the linear skewing function utilized, and various techniques used in number theory. These 
demapping functions which provide the depth measure ak>ng the various storagefretrieval directions, have 
been derived for the 6 and 26 storage/retrieval directions. However, utilizing the principles of the present 
invention, other methods can be derived, for storing and retrieving voxels along more than 26 conflict-free 
storage/retrieval directions. 

In one embodiment of the present invention a transparent parameter which is either transparent or non- 
transparent voxel value is used in order to determine which voxels will quafify for the local competition 
amongst the voxel depth measures. . _4 * 

In another embodiment of the present invention, only those voxel depth measures determined to be 
within a predetermined depth measure range (determined by the observer/user), are allowed to compete in 
the voxel depth competition process. As a result it is possible to selectively "peel away" particular layers of 
35 voxel base images, and determine which voxels are closest on such a basis. 

One aspect of the present invention concerns apparatus for storing and retrieving simultaneously, a 
beam of voxels disposed parallel to a storage/retrieval direction of a n^ matrix anray in 3-D discrete voxel 
space In general, the apparatus comprises storage nneans for storing voxels into memory storage space by 
mapping the voxels along one of the plurality storage/retrieval directions, into a plurality of an independently 
40 accessible memory storage modules. Each memory module is indexed by a memory module index and has 
internal memory cell indices. The mapping is carried out in accordance wltti linear skewing function 
expressed in terms of x. y and z coordinate directions and integer n. A retrieval means is also provided for 
retrieving the voxels along one of tiie plurality of storage/retrieval directions In 3-D voxel space, by 
demapping one or more of the voxels from the memory storage space into tiie 3-D voxel space using 
45 Spatial parameters and integer n. 

In ttie preferred embodiment, ttie storage/retrieval apparatus of the present invention compnses a 
plurality of n modules which are indexed with memory module indices k=0. 1. ...n-1. Each k-th module 
includes the k-th memory module, a k-th local address unit for independently accessing a k-th memory 
module, and each memory cell of each k-th memory module is specified by internal module indices i and j. 
50 Each kHh local address unit Includes a local depth computation unit for computing a voxel deptti measure 
of each voxel residing in the k-th memory module. The voxel depth measure of each voxel is representative 
of the depth of the voxel in tiie n^ matrix array, and is detemnined on the basis of one or more on the x. y. 
and z coordinate values of the voxel, the integer co-efficient a, b. and c of the linear skewing function, and 
the module index k of the memory module into which the voxel has been mapped. 
55 The k>cal addressing unit further includes means for determining for each voxel, the memory cell 
indices i and i. on the basis of tiie voxel depth measure and the x. y and z coordinate values corresponding 
to each voxel. The local addressing unit also includes a memoiy cell selection means for selecting 
cell indices I and j. from global memory cell Indices i and i and the voxel depth measures along with 
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storage/retrieval directions. ^ ■ •, » 

The apparatus ot the present invention also includes a central control unit having a addressing umt for 
generating module indices k" and global memory cell indices i and j . . , ^ 

The central addressing unit further includes means for carrying out the mapping of the voxels along the 
beams parallel to the storage/retrieval direction, into the plurality of an independently accessible memory 
modules. Notably, the mapping is carried out in accordance with the Inear skewing function as discussed 
hereinbefore. Also provided is a means for transferring the global memory cell indices i and j from the 
central addressing unit to the plurality and local address units. In addition a means is provided for 
transferring the beams of voxels into or from the n menrtory modules. 

In the preferred embodiment each k-th module comprises a voxel or multiple write processing unit, a k- 
th voxel value storage register, a k-th voxel depth measure register, a k-th voxel depth competition unit and 
a writing means. The k-th value stored register is for storing the k-lh voxel value of the k-th voxel atong the 
beam parallel to the storage/retrieval direction. The k-th voxel depth measure register is for stonng the voxel 
depth measure of the k-th voxel atong the beam parallel to the storage/retrieval direction. Providing means 
for locally comparing the voxel depth measure of the k-th voxel of ad of the voxels of the beam which is 
parallel to the storage/retrieval directton. and for detennining which voxel depth measure has a predeter- 
mined extreme measure, and in response thereto, generating a voxel pass control signal corresponding to 
the voxel depth measure hawing the predetermined extreme value. The writing means on the other hand, is 
(or writing into a two dimenstonal frame buffer, the voxel value stored in a k-th voxel value storage register 
M which corresponds to a voxel pass control si^. 

The apparatus also includes means tor provWIng transparency control parameters to the voxel multiple 
write processing units. These transparency control parameters are received by a k-th voxel transparency 
control unit, so as to control the kjcal comparison of voxel depth measures within the k-th voxel depth 
competition unit Notably, the voxel depth competitton control is carried out on the basis of the voxel value 
2S in the k-th voxel value register in the k-th transparency control parameter. 

And yet another embodiment <rf the present invention, the storagefretrieval apparatus is indudad within 
the 3-0 voxel based graphic system including: a 3-D geometry processor, a 3-D frame buffer processor, 
and a 3-D viewing processor. The voxel based graphic system comprises a memory and addressing 
system tor storing and retrieving an 3-0 voxel based image along a plurality of storage/retrieval directions 
and 3-0 discrete voxel space. The memory and addressing system includes a plurality of modules, wherein 
each module has a memory module and a local address unit tor independently addressing the memory 
module. A central addressing unit is also provided for simultaneously transmitting to each module, 
broadcasted addresses and control parameters. In addition a transfer means is provided for transferring 
broadcast addresses and control parameters to each module, and in addition a means is also provided for 
transfentng the voxel based bnages Into or from the n memory modules. 

The voxel-based workstation of the present invention is capable of providing a full range of inherent 3-D 
interactive operations In a simple, general work-bench sel-up. since the workstatton operates in both 
discrete 3-0 voxel space and 3-D geometry space, and provides ways in whteh to interact two spaces. 
Accor^gly. the workslaUon can be used with inherent 3-0 user interactive devices, techniques and 
electronic tools, whfeh support direct and natural interaction, generation and editing of 3-0 continuous 
geometttaal models. 3-D discrete voxel Images, and their transformations. Such a 3-D voxel-based 
workstation is appropriate for many 3-D applications such as medical imaging. 3-D compirter-aided design. 
3-D animation and simulatkio {e.g. flight simulation). 3-D image processing and pattern recognition, and 
general 3-0 graphics imeraction. 

A principal advantage of the present inventton is that a 3-0 matrix array of voxels in V space, can be 
mapped into a plurality of independently accessible memory modules with conflict-free access to beams of 
voxels along a plurality of storage/retrieval directions in 2? space. 

Another advantage of the present invention is that beams of voxels along a plurality of conflict-free 
storing/tetrievbig (and viewing) directions in 3-D discrete voxel space 2?. can be stored in and retrieved 
from 3-0 computer storage space, requiring specificatton of the voxel beam in P space solely in terms of 
spatial parameters having a physical interpretation to the user/bbsewer. 

Another advantage of the present invention is that a voxel beam can be mapped from a partcuiar 
confUet-frae storage/retrieval direction in 2? space, into a plurality of independently accessible mwrnt)/ 
modutos in skewed computer storage space according to a linear skewing function which provides conflict- 
SB free access to beams of voxels atong a pluraHty of storage/retrieval directtons. 

A further advantage of the presort inventton is that conlUct-lires access to a voxel beam in the Cubic 
Frame Buffer is provided along storage/retrieval directions parallel to the x. y and z principal aafflS thereof, 
wherein specificatton of the voxel beam, is solely in tenns of spatial parameters, such as the physical 
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coordinates of th© voxels and their physical depth Into the 3-0 matrix array model of the Cubic Frame 
Buffer. 

A further advantage off the present Invention is that conflict-free access to the Cubic Frame Buffer is 
provided along twentysix symmetrically disposed storage/retrieval (e.g. viewing) directions. 
5 Another advantage of the present invention is that the non-transparent voxel in a beam along a specified 
storage/retrieval direction which is closest to an observer along that direction, can be determined. In this 
determination, computed voxel depth measures for voxels along the beam, compete in a process to 
determine the closest voxel. Furthermore, the method can be used to determine which voxel is the closest 
to the observer and which satisfies observer-selected transparency control parameters, such as for 
10 example, color, translucency ratios, texture, transparency and opaqueness. 

A further advantage of the present invention is that depth and interval control parameters are used to 
restrict the voxel depth competition process amongst those voxels falling with ranges specified by such 
depth and/or inten/al control parameter. 

Advantageously, the method of the present invention can be utilized in order to provide a method of 
18 projecting (from 3-0 voxel-based images). 2-0 voxel images along a specified viewing direction. 

An even further advantage is that methods and apparatus utilizing the present invention can be used for 
storing, accessing, and/or processing general-purpose 3-D data in a computer system. 

Other and further advantages of the present invention will be explained hereinafter, and will be more 
particularly delineated in the appended claims, and other advantages of the present invention will be 
20 apparent to one with ordinary skill in the art to which the present invention pertains. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 

For a full understanding of the obiects and advantages of the present invention, reference Is made to 
the following detailed description of the preferred embodiment which is to be taken in connection with the 
accompanying drawings, wherein: 

Fig. 1 is a block diagram of a 3-D voxel-based graphics workstation In which the methods and 
30 apparatus of the present inventk>n are defMcted: 

Rg. 2 is a schematic diagram of a cellular matrix array model of a 3-D Cubic Frame Buffer of the 
voxel-based graphics system of Rg. 1; 

Rg. 3 is a schematic diagram illustrating 2-D and 3-D scan-conversion processes and their 
relationship with discrete 2-D pixel-Image space and discrete 3-D voxel-image space, respectively; 
35 Rg. 4A is a 2-D graphical representation of the "skewed" organization of a 3-D cellular mathx anray 

model of a 4x4x4 Cubic Frame Buffer, showing all voxels on the same diagonal section (i.e. in the same 3-0 
discrete voxel subspace Zj ) modulo 4 being grouped in the same memory module, and also three voxel 
beams in the -x. +2, and +y orthographic directions laid out separately with their respective distribution of 
voxels among the modules, with the internal memory cell addresses (i.e. Indices) shown in parentheses; 
40 Fig. 4B is a graphical representation of the (k=0)-th memory module for the 4x4x4 Cubic Frame 

Buffer of Rg. 4A, showing the diagonal sections of each memory module enctosed In solid line boundaries; 

Rg. 4C is a graphical representation of the (k = l)-th memory module for the 4x4x4 Cubfc Frame 
Buffer of Rg. 4A, showing the diagonal sections of each menwry module enclosed in solid line boundaries; 
Rg. 40 is a graphical representation of the (k = 2)-th memory module for the 4x4x4 Cubic Frame 
45 Buffer of Rg. 4A, showing the diagonal sections of each memory nruxlule enclosed in solid line boundaries; 

Rg. 4E is a graphical representation of the (k=3y-th memory module for the 4x4x4 Cubic Frame 
Buffer of Rg. 4A, showing the diagonal sections of each memory module enclosed in solkl line boundaries; 

Rg. 5A is a graphical representation of a cellular matrix array model of a 3-D Cubic Frame Buffer, 
having six orthographic storage/retrieval (and viewing) directions; 
so Rg. 5B is a graphical representation of a cellular matrix array model of a 3-D Cubic Frame Buffer, 

having twenty-sbc storage/retrieval (and viewing) directions; 

Fig. 6A is a schematic representation of the process of mapping voxel physical coordinates x. y and 
z in 3-D discrete voxel space. Into memory module address indices k and memory cell address indices I 
and j. according to the principles of the present invention; 
59 Rg. 6B is a schematic representation of a process of de-mapping the physical x, y and z coordinates 

of a voxel mapped into the k-th memory module, into a voxel depth measure, d. along a particular 
storage/retrieval (i.e. viewing) dlrectton. and subsequently into the memory address indices k. I, j. according 
to the principles of the present invention: 

8 
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Fig. 7 is a block system diagram of the 3-0 voxel-based workstation constructed In accordance with 

P^^P-J-' -^-^^^^^ 3.0 .<,.e.based workstation o, the ^sent ;nve««ojs.2^^^^ 
Drocessors a central addressing unit, and a plurality ot n modules each having a local addressing unit an 
s toidependently accessible memory module, and a voxel multiple write processing -nd 
RQ 9 is a block schematic diagram o« a k-th module of the apparatus for stonng. accessing and 
processing voxel-based data according to the principles of the present InvenHon: 

10 is a block schematic diagram of the memory and addressing system of the present invention 
Illustrated in Rg. 9. showing the central addressing unit and several modules; 
,0 Rg. Ills a schematic representation of a portion of the voxel muHlpte wnte 

the conjuration of the k-th voxel depth computation unH of the tocal 'j'^^^'^J^^ 
mLiS «Lage register, the k-th local voxel depth competition unit, and the voxel depth bus of the present 

^""^g. 12 is a graphical representation of the Intermodule voxel value moving (l.e. shifting) mechanism 
« ofthe gesent ;-n«o^ -d^^^^ ^ ^^^^ „^ ^ ,ocal addr^g un^ 

and k-«? memory module of the apparatus for storing, accessing and processing voxeH»ased data 
according to principles of the present invention. 
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Referrina to Figs 6A. 68. and 5A and 5B in particular, the method of storing and retrieving beams of 
v a pluSity oi storage/retrieval directions in 3-D discrete voxel space, according to the present 

'"^T^s onSS, ?h?«se Of Six storage^etrieva. direCons ^^1 x. ^J^^^'^ 
axes^B beconsidered. in Fig. 5A. a 3-0 matrix array model of the Cubic Frame Buffer In 3^) ««• fP«« 
^te fflustSTwhe^in tfi Six orthographic storage/retrieval directions in ? space are shown by 
directional arrows labelled si. t2 and t3 for the tx. ty and « principal axes, respectively. 

In ? voxel space, each voxel has physical coordinate values x. y and z ''"^ ' 
directi^ respectively, which taken together form a 3-0 matrix array. Passing through tfie 3-0 "«Ji>. «na^ 
aSgs S^els referred to hereinafter as "beams-, can also be called "rays" as well de^l.^ on »« 
^tf ?mptementation involved. ^ tne case where the Cubic Frame BuHer is a memory cu^ of n3 voxe s 
Sraddress indices of each voxel V{x. y. z) in a beam of n voxels B,. are represented as k. .. i m 3 D 
comouter storage space N'. which will be further specified below. . • «a»,iv 

Tn order to enable simultaneous (i.e. conflict-free) access to a full beam <>' " v««rts in *e 3^ m^.x 
arrav in Z3 space regardless of the storage/reUieval direction, the large number of memory banks Q.e. 
Sisl^SS tl^Cubic Frame Buffer 2. must be accessed independently in parallel. It is just th^sort 
7^l:^ZlX^U, that is required in order to .^^^^^^'^^^^ 
technique involving processing a full beam of voxels at a time, .n order to seteet^voxel atong each beam 
which is closest to an observer viewing along a viewing direction parallel to the beams. 

1\,e efficiency of such a conflict-free memory access method for "voxel beam." processmg 

may^iSTSom the .act that in one memory cycle, a full beam of r, voxels could ^T^^^^^^^ 
!e element) from each of the n memory banks provided in the memory storage and accessing system 
Slower it s aTon Sal problem to store the voxel elements of a 3-0 matrix array such that all beams of 
Se.i alonSL x. y and z principal axes) can be stored and retrieved in one memory access cycle. 

i.e. as beams of n voxels, "without conflict". _f-„~,a 

^rxel beam storage and retrieval method of the present Invention relies on a memonr storage 
schemJ ^^t rn J^he voxels of a 3-0 matrix array (Cubic Frame Buffer) in Z3 space^into a certan num^^ 
rrrJoHTmod^les ,i.e. banks) and provides tor the conflict-free access to 

to a plurality of storage/retrievd directions. Such a memory storage is <^^^^J^ JlfSSS 
Morage scheme utilized. Notably, the practicality of skewir,g schemes, .n general, depends upon whether 
fheycan be described by a small amount of tabular information, or by a simple formula, 
"^''k^e method of Jpresent invendon. a gene«l class of '^[^"^J^"^'^^^ 
are used in mapping a 3-D matrix array of voxels (e.g. a 3-0 voxeW»ased irnage) into a pluralrty ofn 
Li^moduirindexed k=0. 1. ... n-l. Each k-th memory module corrtatos n» memory cells, with each 
memory ceU being spedfted by internal memory module mdices i and j. 
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The general form of the linear skewing function S(x, z) is given by 

S (X, y, 2) « k = (ax + by + cz) mod n 
i « X j * y 

where x, y and z are the coordinate values of each voxel in 3-0 discrete voxel space Z^, 0 ^ k, x, y. z £ n-1, 
where a, b and c are integer coefficients, and where n is a prime relative to the integer coefficients a. b and 
c and to linear combinations thereof. e.g. a - b. b c, a b -t- c, etc. 

In general, the linear skewing function S maps physical coordinate values x, y and z represented in 3-0 
voxel space Z^. into memory module storage Indices k, i. j represented in 3-D skewed memory storage 
space N^. This mapping process is illustrated in Rg. 6A. in particular. 

In order that conflict-free access of voxel beams will be provided parallel to particular storage/retrieval, 
or viewing directions into the 3-0 matrix array (i.e. Cubic Frame Buffer), it is essential to formulate the 
condlttons on a. b. c and n such that linear skewing scheme S(x, y, z) is conflict-free on. e.g. the 
orthographic tx. ±y, and tz axes of the 3-D matrix array in space. As will be explained more fully 
hereinafter, conflict-free access to specified storage/retrieval viewing directions can be provided for a proper 
choice of a, b. c and n. which can be achieved by ensuring that no two voxels along a beam parallel to a 
storage/retrieval direction, are mapped into the same Independently accessible memory storage nnodule K. 

Refening now to Rg. 5A. the six storage and retrieval directk)ns of the method of the present invention 
are specified as foltows. "Rows" are referred to hereinafter as those storage/retrieval and viewing directions 
that are parallel to the x principal axis. "Columns" are referred to hereinafter as storage/retrieval directions 
that are parallel to the y principal axis, whereas, "axles" are referred to as storage viewing directions that 
are parallel to the z prindpal axis. 

By utilizing a 3-D memory device comprising n independently addressable memory modules having n^ 
memory cells each, it is clear then, that confllct-firee access to t>eams of voxels parallel to any one of the 
six storage/retrieval directions specified in Fig. 5A, will be provided by guaranteeing that no two voxels in a 
beam parallel to one of the storage/retrieval directk^ns reside in that same k-th memory module. 

To more fully appreciate the method of the present invention, it will be helpful to briefly discuss several 
basic concepts in number theory, for example, the concept of the greatest common divisor, and the concept 
of "relatively prime", are in order. 

Every integer which cfivides all the integers a. b I is said to be a "common divisor** of them. The 

largest of these comnton divisors is said to be their "greatest common divisor" (i.e. gcd) and is denoted in 
the mathematics of number theory, by the symbol (a, b, .... l)=l. and hereinafter by a somewhat more 
redundant notatkm: god (a, b. ... 1). In view of the finiteness of the number of common divisors, the 
existence of the "greatest common divisor" is evident In this regard, a few definitions are in order. If gcd 
(a, b, ... I) ' 1, then a, b, ... I are said to b9 "relatively prime", which signifies that there is no divisor, 
except for 1, which is common to a, b I. 

If each of the numbers a, b, ... I is relatively prime to any other of them, then a, b I are said to be 

"pairwise prime". It is thus evident that pairwise prime numbers are also relatively prime, and in the case of 
two numbers, the concepts of "pairwise prime" and "relatively prime" coincide. 

In the prefenred embodiment hereof, the linear skewing scheme S(x. y, z) is "proper' in that all memory 
modules are used in the mapping operation. The "proper" condition of the linear skewing scheme S(x. y. z) 
is satisfied if and only if Integer coefficients a. b and c are "relatively prime" with respect to n. This 
"relatively prime" condition for proper linear skewing sciiemes can be expressed altematively as follows: 
The greatest comnnon divisor 0-e. "gcd") of integer coefficients a. b. c and n must be equal to one. i.e. god 
(a.b.c,n) = 1, 

According to such a linear skewing scheme, the voxels of a beam along a row of the 3-0 matrix array 
(e.g. parallel to the -t>x principal axis, with y and z fixed constant) are placed or "ordered" in successive 
memory modules which are "a" apart mod n, and are termed a-ordered voxels. Similariy. the voxels of a 
beam ak)ng a column of the 3-0 matrix array (e.g. parallel to the -i-y principal axis, with x and z fixed 
constant are placed in successive memory modules which are b apart mod n, termed b-onjered voxels. 
Also, those voxels of a beam akxig an axle of the 3-D matrix array (e.g. parallel to the -•'Z principal axis, 
with X. y fixed constant) are placed in successive memory modules c apart mod n, and are termed c- 
ordered. 

For a proper choice of a, b. c and n. the linear skewing functkMi k « (ax «*• by -^ez) mod n provides 
cortflict-free access to all rows, columns and axles of the 3-0 matrix array. Having specified the a, b and c 
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ordering of the voxels along the tx. ty and tz principal axis of the 3-D matrix array in Fig 13A. and the 
condition gcd (a. b. c. n) = 1. these specifications suggest the selection of a, b. c and n as follows: 

CoJi^ntli. n can be any number. e.g. preferably any power of 2 as this simplifies memory construction, 
addressing, and/or the modulo computation in the linear skewing function S(x. y.z). Notably, other selections 

of a. b, c and n can satisfy the condition gcd (a,b.c,n) = 1. . „ « ^. 

Thus, the linear skewing function of k = (x + y * 2) mod n. performs a mapping of a 3-0 matrix array 
in r space, into a plurality of n memory modules of memory cells in 3-0 computer storage space N . 
whereby six conflict-free storage/retrieval directions along the principal axes are provided. 

After having mapped (i.e. stored) beams of voxels (characterized by x. y and z coordinates) from Z 
voxel space into 3-0 computer storage space (characterized by memory storage cell address indices k. 
i i) using the above-derived proper linear skewing function, a user/observer typically has a natural need to 
retrieve the stored beam of voxels from 3-D computer storage space N^. by demapping the beam of voxels 
in 3-D voxel space 2^ for processing and/or viewing. However, as discussed hereinbefore, this poses a 
problem regarding voxel beam retrieval, namely, that individual voxels of a beam once accurately 
represented in discrete voxel space 2? have been scrambled and stored away in n different mernory 
modules having nxn memory cells with address indices i. i that differ from the physical x. y and z 
coordinates of the voxels. In short, the user/viewer is incapable of specifying the memory storage addresses 
corresponding to a particular voxel beam which is to be retrieved along a particular conflict-free 
20 Storage/retrieval direction. . . 

Reference is now made to Fig. 6B in particular, wherein a "demapping" process is illus rated for 
retrieving a voxel beam mapped into skewed memory storage space H>. according to ttie prinaples of the 

'''^TgeSrapplicants have discovered an effective way in which to overcome the aboveKlescrlbed 
voxel-beam retrieval problem. According to the present inventton. a whole beam of voxels mapped into 
computer storage space according to a particular linear skewing function, can be retrieved from skewed 
storage space by a process demapping the voxel beam back into 3-0 voxel space V. Advantageously, 
such a demapping process only requires specifying the beam of voxels in discrete voxel space solely m 
terms of spatial parameters having a physical interpretation to the user/observer. 

in the preferred embodiment, the spatial parameters that have been selected lor voxel-tjeam retrieval 
specification are (i) the physical coordinate direction of ttie voxel beam, and (il) the physhal depth of the 
beam into the 3-D matrix array model of the Cubic Frame Buffer. These spatial parameters "^f 
that they have a physical interpretation, namely they specify the direction of the voxel beam, and also the 
depth of the voxels in the beam extending into the 3-0 matrix array along the specified storage/retneval 

^ ^*'^/te?rrdi?g*t^"method of the present invention, the demapping process involves specifying along 
which of the six storage/retrieval directions, the voxel beam (to be retrieved) extends into the 3-D nwtt.x 
array in 2? space. For ttie case of retrieving the voxel beam along ttie *x storage/retrieval direction, ttie 
user/obsenrer needs only to specify a constant value of y and z. e.g. y.. z,. respectively and compute 
deptti measure of each voxel (referred to hereinafter as -voxel deptti measure d") on the basis of spatia 
paTameters. such as x. y and z. and ttie memory module index k » 0. 1. n-1. Notably h<^«-.«;- 
which is. for example, five (5) units deep into the 3-D matrix array in voxel space Z». ,s not in tlw 
memory module k=5. but has been mapped into some ottier memory module according to ttia linear 

« *^l5r!rc!!i'pute*tt« voxel depth measure values d of the beam parallel to the *x coordinate 
direction, a demapping function is derived from ttie linear skewing function S(x. y. z) = (x + y + c) mod n 
by manipulating the same using known -modular arithmetic" principles and operations. For a d'scussion 0 
the invoked principles and operations of modular arithmetic, reference is ma<te to Bemwrts ^ 
Theory by I.M Vinogradov, translated from ttie Fifth Revised Edition (1949) by Saul Kravetz. published by 

so oS^Publfeations in 1954. The result of ttie manipulation of ttie linear skewing function is a demapping 
function which provides the voxel depth measure d. atong Oie +x storage/retrieval direction as folkiws: 

a (k • (y + 2)) mod n. ^ , a < «.i 

Thus by evaluating in parallel, ttie voxel deptti function d« for all values of memory module k = 0. 1. ... fM. 
ttie memory storage addresses, or indices, for the entire beam to be retrieved (i.e. accessed) along the 
as specified storagefretrieval direction, are given by i = dx. i = y for k » 0. 1. 2 ... n-1. 

When retrteving the voxel beam along ttie -x storage*etrieval direction, the corresponding voxel deptti 
measure parallel to ttiat direction d« can be computed by 
d, ■» n -d»,- 1. 
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In a similar fashion as described above, demapping functions for providing voxel depth measures of 
voxels of beams along the ty and tz coordinate directions can be derived in the form as follows: 
d^y « (k-(x ♦ 2)) modn 
5 d.y = n - d*y + 1 

dz = (k - (x + y)) mod n 
d.2 = n-d»^-l. 

By evaluating in parallel, the voxel depth functions d^y or d^z, for all values of memory modules indices 
k a 0. 1, ... n-1. the memory storage addresses for the beam to be retrieved along the specified 

10 storage/retrieval direction, are determinable in a manner described at>oye. 

Refemng now to Fig. 6A, 6B. and Fig. 58 in particular, the case of twenty-six conflict-free 
storage/retrieval directions of a 3-0 matrix array in 3-0 voxel space Z^. will be considered. 

In Fig. SB. a 3-0 matrix array model of the Cubic Frame Buffer in 3-0 voxel space is illustrated, 
wherein the twenty-six storage/retrieval directions in space are shown by directional arrows labelled t1 

75 for the rows parallel to the tx principal axes. t2 for the columns parallel to the t principal axes, and t3 for 
the axles parallel to the ±z principal axes. The other additional twenty storage/retrieval directions are 
specified along the minor-diagonals represented by directional an^ows ±4. ±5. £6. the minor-anti-diagonals 
represented by directional arrows £7. ±8. ^9. major-diagonals represented by directional arrows tio. £12. 
and major anti-diagonals represented t)y directional arrows £11. £13. 

20 As in the case of six conflict-free storage/retrieval directions, the voxel beam storage and retrieval 
method of the present invention relies on a memory storage scheme that maps the voxels of a 3-D matrix 
array (Cubic Frame Buffer) in 3-0 voxel space Z^. into n independently accessible memory storage 
modules, which provides for conflict-free access to beams of voxels parallel to any one of twenty-six 
storage/retrieval directions. Again, the storage mapping process utilizes a linear skewing functton of the 

25 following general form: 
k = (ax + by + cz) mod n 
where 0 £ k, x. y. z ^ n-1. 

To ensure the conflict-free access to rows, columns and axles, the condition gcd (a. b. c. n) = 1 must 
be maintained. For a proper choice of a. b. c and n. the above-skewing function will map the voxels of a 3-D 

30 matrix into n independently accessible memory modules, so as to provide conflict-free access to ail rows, 
columns, axles, major and minor diagonals and anti-diagonals. 

In order to provkje a storage scheme such that a beam of voxels can be retrieved conflict-free from 
memory storage space, along any one of the twenty-six storage/retrieval directions in V space specified 
above, it has been discovered that the following voxel ordering conditions must be satisfied, namely: 

35 The voxels of a row are placed (i.e. mapped into, or "ordered") in successive memory modules which 
are "a apart modulo n". termed a-ordered. The voxels of a column are placed in successive memory 
modules "b apart modulo n", termed b-ordered. The voxels of an axle are placed In successive memory 
modules "c apart modulo n", termed c-ordered. The voxels of major diagonals are placed in successive 
memory modules "(a + b + c) apart modulo n", termed (a + b + c)-ordered, and "a ♦ b - c apart 

40 modulo n", termed (a b - c)-ordered. The voxels of a maRMr anti^iagonal are placed In successive 
memory modules "(a + b - c) apart modulo n", termed (a + b - c)-ordered. and "a - b* c apart modulo 
n". termed (a - b c)-ordered. The voxels along minor diagonals and anti-diagonals are placed in 
successive memory modules "(a + b) apart modulo n", termed (a + b)-ordered. "(a - b) apart modulo n" 
termed (a • b)-ordered in the x«y faces, are (a -*> c)-ordered. and (a - c)-ordered in the y-x faces, and the 

45 voxels are temied (b * c)-ordered and (b • c)-ordered in the y-x faces. The voxels along the positive and 
negative directions of the same beam have the same "order". 

In view of the above voxel-ordering conditions that must be satisfied for the case of twenty-six conflict- 
free storage/retrieval directions, the selection of integer coefficients a. b. c and n are restricted. Such 
restrictions can be formulated in terms of Bask: Concepts of Divisibility Theory, and in particular, in terms of 

50 the concept of the "Greatest Common Divisor" (i.e.. gcd). the theoretical foundations of which are discussed 
in detail in Chapter 1 . Section 2 of Elements of Number Theory by I.M. Vinogradov, supra. 

The restrictions Imposed on a, b, c and n for the case of the above specified twenty-six storage/retrieval 
directions of tiie 3-D matrix anray in voxel space Z?, are given as folbws. 

55 
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a,*b a^c bfc a + b7»c a7»b + c b^a + c)(l) 
god (a,n) =1 9^ M =1 ^od (c,n) - 1 ) (2) 



qcd (a + b -w, n) = 1 god (a + b - c, n) - l 
god (a - b -c, n) » 1 god <a - b + c, n) - l 
god (a + b, n) = 1 god (a - b, n) = 1 



) 

) (3) 
) 

) (4) 



god (a ♦ c, n) - 1 



god (a - c, n) » 1 > (5) 



god (b + c, n) - 1 god (b - c, n) - 1 



) (6) 



For example, a =5. b =2 and c = 1 satisfy equations (1). and thus the general slowing function becomes 

k = (5x + 2y + 2) mod n. ,« u. u 

" fstotably integer n is selected so as to satisfy the greatest common divisor equations (2) through (B). 
Accordingly, to realize a 5123 voxel Cubic Frame Buffer 2. n can be selected, for example, from integers 
n-S09 or n«5l7. which are both "prime relative to" integers c = 1 through integer (a+b+c)=8. as 
required by the relatively prime conditions formulated in equations (2) through (6). 

Having now specified one set of integer coefficient values lor a. b. c. and n in the generalized hnear 
skewing function S(x. y. z). the voxels in the 3-D matrix array in voxel space can be mapped into 
memory storage space as discussed hereinbefore. Once mapped into storage space, the voxels or 
beam thereof along one ol the twenty-six conflict-free storage/retrieval directions m space, can be 
retrieved by a demapping process involving demapping of the voxels from storage space into voxel 
space using spatial parameters and integer n. The details of this demapping operation in the voxel beam 
^ retrieval process, will now be described below in detail for the case of twenty-six confUct-free 
storage/retrieval directions. 

As with the case of six conflict-free storage/retrieval directions, retrieval of one or more voxels of a 
beam in 3-0 voxel space Z-'. requires, according to the principles of the present invention, determination of 
the memory storage addresses fin 1^ storage space) of the voxels of the beam to be retrieved parallel to 
one of the twenty-six conflict-free storage/retrieval directions. 

However, since the voxels stored away in memory storage space are "scrambled as a result of the 
skewed storage mapping, determination of the memory storage addresses of the voxel beam is not 
straightforward, and requires that a different demapping operation be carried out for retrieval along each ot 
the twenty-six storage/retrieval directions. Thus, specification of those demapping operations for thirteen 
viewing directions Is provided betow for positively specified storage/retrieval directions. Notably, the 
demapping operations for the negatively specified "antiparaller storage/retrieval directions are denvable in 
a stratghtfonward manner. Fdr purposes <^ illustration, integer n representative of the number of indepen- 
dently accessible memory modules, is taken to be 517. for the linear skewing funcUon: 
S (x. y. 2) a k a (5x + 2y + z) mod 517 

*^*^ln~order to access a beam of voxels along "rows" in the 3-0 matrix array (i.e. parallel to the x-axis). 
storagefretrieval direction No. + 1 illustrated in Fig. 13B is uBIized. As illustrated in Fig. 13B. each row in the 
3-0 matrix array is defined by coordinate y « y, and coordinate z = z,. Along this storage/retneval 
direction, internal memory cell address j = y « y^, and each memory module ad*ess k along the row 
" specffies a different i »x. which is the voxel depth measure d.. i.e. distance from the observer along the + x 
axis. 

To derive a mathematical expresston for d., which is a demapping function, the Bnear skewingji^ 
mapping) function k => (5x ♦ 2y + z) mod 517 is manipulated using modular arithmetic operattons. (which 
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are different from conventional algebra), to obtain the following congruence: 
5x = k-(2y + 2) mod 517. 

In order to obtain x alone on the left side of the above expression, which represents the voxel depth 
measure (i*n, resort must be made to the fundamental properties of fractions, and in particular the theory of 
5 Continued Fractions, to which reference is made in Chapter 1. Section 4 of Elements of Number Theory . 
supra . 

The ratio of n/a, e.g. 517/5, can be developed into a continued fraction n/a. 
Considering the last two convergents. 



10 



m 



n/a, 



15 

by the properties of congruences in one unknown (discussed in Chapter IV, Sectton 2, of Elements of 
Number Theory , supra) , the following congruence of the general form 
A X s B mod n 
can be expressed as 
20 X » «-l)'"PMB)modn 

which, for convenience of notation, will be represented as 
X s <T (5 ) • B) mod n 

Thus, using the above moduli relationships, the congruence 5x k - (2y -i- z) mod 517 can be expressed 
as 

25 X = 207 (k - 2y - 2) mod 517 = d«. 

517 

Where T ( -g- ) = 207. 

Accordingly, using dx. j = y = yo and z = Zo for k = 0. 1. ... n-1. all of the addresses k. i, j of the voxels 
along the >x direction are specified, and thus the full beam of voxels can be retrieved along the 
30 storage/retrieval direction No. -i-l (illustrated in Rg. SB) in terms of spatial parameters. i.e. yo and Zq, and 
integer k. 

To detemnine the voxel depth value d.^ along the -x axis, d.^ is simply computed by the following 
formula: 

d-x — n - d*x • 1- 

35 For the general case of voxel beam retrieval along storage/retrieval direction No. +1, where the 
skewing function k « (ax by + cz) mod n is employed in storage mapping operation, the derived 
congruence 

ax = (k - by - cz) mod n 

can be manipulated as hereinbefore described, to provide the voxel depth measure d«x as 
^ d*x = X » (Tg ) (k - by ' cz)) mod n. 

In order to access a t>eam of voxels along "columns** in the 3-D matrix array (i.e. parallel to the y-axis 
thereof), storage/retrieval direction No. *2 illustrated in Rg. SB is utilized. As illustrated in Rg. SB. each 
column in the 3-0 matrix an-ay is defined by coordinate x=Xo and coordinate z^z©. 

For the general case of voxel beam retrieval along storage/retrieval direction No. +2. where the 
^ skewing function k = (ax + by + cz) mod n is emptoyed in storage mapping, the derived congruence 
by » k - (ax •*> cy) mod n 

can be manipulated as hereinbefore described, to provide the voxel depth measure d«y as 
d*y = y = T (5 ) (k - ax - cz) mod n. 

In order to access a beam of voxels atong "axles" in the 3-D matrix anray (i.e. parallel to the z-axis). 
5^ storage/retrieval direction No. +3 as illustrated in Fig. SB is utilized. As illustrated in Fig. 5B, each axle in 
the 3-D matrix array Is defined by coordinate x^Xo. and coordinate y =yo. 

For the general case of voxel beam retrieval along storage/retrieval direction No. -(-3. where the 
skewing function k » (ax by cz) mod n is employed in the storage mapping operation, the derived 
congruence 
55 cz = (k - (ax + by)) mod n 

can be manipulated as hereinbefore described, to provide the voxel depth measure d.x as 
d*, » 2 « (Tg )(k - ax - by)) mod n 

Accessing a beam of voxels along *^e minor diagonal" in the 3-D matrix array specified by 



14 



BNSDOCID; <EP 0368425A2 I > 



EP 0 368 42S A2 



storage/retrieval direction No. +4. is now considered. As iUustrated in Fig. 5B. each minor diagonal in the 3- 
D matrix array is defined by equation x-y = constant and coordinate z = z„. The fact that x-y = constant 
provides two options tor determining the depth along the beam. According to the first option, along this 
storage/retrieval direction, each memory module address k specifies a different coordinate y, which is 

s proportional to the distance dy from the obsewer along the minor diagonal. Alternatively, each memory 
module acWress k also spedfies a different coordinate x, which is proportional to the distance dK from the 
observer along the minor diagonal. 

The case where each k along the minor diagonal specifies a different y coordinate, will be considered. 
Using a procedure similar to that employed hereinbefore, the depth measure d*y Is obtained as 

to d*y = y = (74 (k • 5(x-y) • z)) mod 517 

517 

forn = 517. a = 5. b»2.c = 1.andT(-Tp- ) =74. 

For the gaieral case of voxel beam retrieval along storage/retrieval direction No. +4. where the skewing 
function k = (ax ♦ by + cz) mod n is employed in the storage mapping operation, the derived congruence 
;s k = (ax - ay + (a ♦ b)y + cz) mod n 
can be manipulated to provide 
(a ♦ b) y = (k • (a (x-y) + cz)) mod n 

from which the voxel depth measure dy can be derived using modular arithmetic: 

20 d*y = y=s(T( " )(k-a(x-y)-cz))modn 

Alternatively, for tt^second option discussed above, d^ Is provided by: 

d.x » X = (T ( ) (k + b (x-y) • cz)) mod n. 

Accessing a beam of voxels along the minor diagonal specified by storage/retrieval direction No. +5 in 
the 3-D matrix array, will now be considered. As illustrated in Fig. 5B, each such minor diagonal is defined 
by equation y-z = constant, and coonJinate x=Xo. 

For the general case of voxel beam retrieval along storage/retrieval direction No. +5. where ttie 
generalized linear skewing function is employed in the storage mapping operation, the derived congruence 
k 3 (ax + by • bz + (c + b)z) mod n 
can be manipulated to provide 
(c + b)2 = (k - (ax + b(y-z))) mod n. 

from which the voxel depth measure d«, can be derived using moduli algebra, as fotows: 

d»,«z = (T( ) (k • ax - b{y-z))) mod n 

3S Alternatively, for the second option discussed above, d^y is provided by: 

d*y = y = (T( ) (k-ax c(y - z))) mod n ^. . ^. « - 

Accessing a beam of voxels along the minor diagonal specified by storage/retrieval direction No. +6 in 
the 3-D matrix ^ay will now be considered. As illustrated in Fig. SB, each such minor diagonal is defined 
by ^uation x-z ^ constant, and coordinate y - yo- 

For the genwal case of voxel beam retrieval atong the storage/direction No. *6, where the generalized 
linear skewing function is employed in the storage mapping operation, the derived consequence 
k = (ax • az + by + (c +a)z) mod n 
can be manipulated to provide 
(c + a)2 = (k - (a (x-z) + by)) mod n. 

from which the voxel depth measire d* j can be derived using modular arithmetic, as follows: 

d*i » z « T (( -57^ ) (k- a(x-2) - by)) mod n 
so Alternatively, for tho second option discussed above, d^x is provided by: 
<*-» = n'(-i7r )(k-by * c(x.z)))modn. 

Accessing a beam, of voxels along the minor anti-diagonal specified by storage/retrieval direction No. 
+7 in the 3-D matrix array, will now be considered. As Illustrated in Fig. 5B. each such minor anti-diagonal 
is defined by the equation x ^ y ^ constant, and coordinate z » z^. 
65 For the general case of voxel beam retrieval along storage/retrieval direction No. +7. vifhere the 
generalized linear skewing function Is employed in the storage mapping operation, the derived congruence 
k a (ax ♦ ay + (b-a) y ♦ cz) mod n 
can be manipulated to provide 
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(b • a) y = (k - (a + y) + cz)) mod n. 

from which the voxel depth measure d^y can be derived using moduli aigebra. as follows: 
n 

d*y = y = (T( -gzi" ) (k - a <x -i-y).- cz)) mod n. 
5 Alternatively, for the second option discussed above, d^x »s provided by: 

n 

d*x - K = (T( -51:5- ) (k - b (x+y) . C2)) mod n 

Accessing a beam of voxels along the minor anti-diagonal specified by storage/retrieval direction No. +8 in 

the 3-D matrix array, will now be considered. As illustrated in Fig. 5B. each such minor anti-diagonal is 
10 defmed by equation y -»> z = constant, and coordinate x s Xq. 

For the general case of voxel beam retrieval along storage/retrieval direction No. +8. where the 

generalized linear skewing function is employed in the storage mapping operation, the derived congruence 

k - (ax + by + bz + (c-b)2))mod n 

can te manipulated to provide 
f 5 (c-b)z = (k-ax-b(y + z)) mod n 

from which the voxel depth measure d*, can be derived using modular algebra, as follows: 

d+2«z=(TY^()c-ax-b(y+2))) mod n 



Alternatively, for the second option discussed above. d«.y is provided by: 

d+y-y»(T^^r5j-(lc-c(y+z)-ax)) mod n 

Assessing a beam of voxels along a minor anti-diagonal specified by storage/retrieval direction No. +9 

in the 3-D matrix an'ay. will now be considered. As illustrated in Rg. 5B. each such minor anti-diagonal is 

defined by the equation x = z - constant, and coordinate y = yo- 
^ For the general case of voxel beam retrieval along storage/retrieval direction No. +9. where the 

generalized linear skewing function is employed in the storage mapping operation, the derived congruence 

ks5(ax+a2 + by + (c-a)2) mod n 

can be manipulated to provide 

(c-a) z = (k-(a(x + 2)+by)) mod n 
^ from which the voxel depth measure d*, can be derived by using moduli algebra, as follows: 

d*a = z « (T( ) (k.a(x+z)-by)) mod n 

Alternatively, for the second option discussed above. d*x is provided by: 

^ d*« = x=(T( - 3^0 ■ ) (k-c(x+z)+by)) mod n 

Accessing a beam of voxels along a major diagonal specified by storage retrieval direction No. +10 in 
the 3-D matrix array, will now be considered. As illustrated in Rg. 5B. each such major diagonal is defined 
by the equations x-y* constant, y-2« constant and x-zs constant, (each one of these equations can be 
derived from the other two). These equals provide three options for determining the "depth measure" along 
the beam. According to the first option, each memory module address along the major diagonal specifies a 
different coordinate x which is proportional to the distance (i-e.. depth measure) d^^ from the observer atong 
the major diagonal. 

Altematively. each memory module k also specifies a different coordinate y which is proportional to the 
distance measure d*y from the observer along the major diagonal. In the third option, each memory module 
^ k specifies a different coordinate z which is proportional to the distance measure d*r from the observer 
atong the major diagonal. 

The case where each k along the major diagonal specifying a different x coordinate wilt be considered. 
Using a procedure similar to that employed hereinbefore, the voxel depth measure d^x >s obtained as: 
d*K=x=»(194(k+(y-z)+3{x-y))) mod 517 
^ For the general case of voxel beam retrieval atong storage/iretrieval directton No. +10. where the 
generalized linear skewing function is employed in the storage mapf»ng operatton. the derived congmence 
k«(ax-ay +(b+a)y-(b+a)z+(c+b+a)z) mod n 
can be manipulated to provide 
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(a ♦ b + c)2 3 (k-(a(x.y) + (b + a)(y2))) mod n 

from which the voxel depth measure d*i can be derived using modular arithmetic, as follows: 
d-^,«2=(T(-Ji— ) (k-a(x-y)-(a+b) (y-z))) mod n 

Alternatively, for the second option discussed above. d*K is provided by: 

to 

dj.^«x«(T(— il— ) (lc-b(y-2)-^(a+b) mod n 

a-i-b+c 

'5 Attematively. for the third option discussed above. 

d+v«y-(T(-Jl-) (K-a(x-2) + (a+c) (y-2))) mod n 
y a+b+c 

20 

Accessing a beam of voxels along a major anti-diagonal specified by storage/retrieval direction No. + 1 1 
in the 3-0 matrix array, will now be considered. As illustrated in Rg. SB. each such major ant-diagonal is 
defined by the equations x+y»conslanl. y^z-constant and x-z»constant (each one of these equations 
^ can be derived from the other two). ^^.^ 

for the general case of voxel beam retrieval along storage/retrieval direction No. +11. where the 
generalized linear skewing function is employed in the storage mapping operation, the denved congruence 
k a (ax ♦ ay ♦ (b-a)y •(b-a)z + (c-b + a)z) mod n 
can be manipulated to provide 

(a-b ♦ c)z a (k-(a(x + y) (b-a)(y ♦ z))) mod n ^ . . 

from whtoh the voxel depth measure d., can be derived using modular arithmetic, as follows: 

d*,=2«(T( IL) (k-a(x+y)-(b-a) (y+z))) mod n 

a-b+c 

35 

Alternatively, for the second option discussed above. d*y is provided by: 

40 d+v»x=i(T( JL) (Ic-b{y+2) + (c-b) (x-z) ) ) mod n 

* a-b+c 

Alternatively, for the Uiird option discussed above. d*y is provide by: 

^ Accessing a beam of voxels along a major anti-diagonal specified by storage/retrieval direction No. ♦ 12 
in the 3-D matrix array, will now be considered. As illustrated in Fig. 5B. each such major anti-diagonal is 
defined by the equations x-y-constant y*zaconstant and x*z«constant. (each one of these equations 
can be derived from the otiier two). ^ 
For the general case of voxel beam retrieval along storage/retrieval direction No. ♦12. where the 

^ generalized linear skewing function is employed in the storage mapping operation, the derived congnience 
k»(ax+az^(c-a)z*(c-a)y+(b-c+a)y) mod n 
can k>e manipulated to provide 
(a^b-c)ya(k-(a(x*2)^(c-aKy*2))>n«ai n 
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from which the voxel depth measure d^y can t>e derived using modular arithmetic, as follows: 

d+^yo(T( 2 — )(k-a(x+z)-(c-a)(y+2))) mod n 

~ -a+b-c 

Alternatively, for the second option discussed above, d^^ is provided by: 

iLv»x«(T(— II— ) (-k-Hb(y+2)-(b+c){x+z))) mod n 
-a+b-c 

Aitematively. for the third option discussed above, d.^ is provided by: 

d.^=gs(T( ^ ) (-k+a(x-y)-(a+b) (y+z) )) mod n 
a+b-c 



Accessing a beam of voxels along a major anti-diagonal specified by storage/retrieval direction No. +13 
in the 3-D matrix anay. will now be considered. As illustrated in Rg. 5B. each such major anti-diagonai Is 
defined by the equations x+y= constant, y-z» constant and x+z*constant, (each one of these equations 
can be derived from the other two). 
25 For the general case of voxel beam retrieval along storage/retrieval direction No. +13. where the 
generalized linear skewing function is employed in the storage mapping operation, the derived congruence 
k = (bx + by + (a-b)x + (a-b)z + (c-a + b)z) mod n 
can be mantpuiated to provide 
(-a + b +c)y = Oc(b(x + y) + (a-b)(x + z))) mod n 
^ from which the voxel depth measure d*z can be derived using modular arithmetic, as follows: 

^j, 3«(Y(_!L-.)(lc-b(x+2)-(a-b)(x+2))) mod n 
^* -a+b+c 

Aitematively. for the second opfion discussed above, d^x is provided by: 

d+v=x-(T( —) (-k+b(y-z)-(b+c) (x+z) ) ) mod n 

40 ^ -a+b+c 

Alternatively, tor the third option discussed above, d»y Is provided by: 



45 



d.z«y=(T( 2— )(k-a(x+z)-(a-c)(y-z))) mod n 

^ -a+b+c 



While only the cases of six and twenty six conflict-free storage/retrieval directions have been herein- 
^ described, it is within the full understanding of the present invention to provide methods of and apparatus 
for storing and retrieving voxels along mae than twenty-six storage/retrieval directions, utilizing the 
principals of the present invention. In particular, more than twenty-six storage/retrieval directions are 
obtainable with the present invention, and such directions can be either "conflict-free" single access, or 
involve multiple but a limited number of conflict-free accesses to the Cubic Frame Buffer in 3-0 discrete 
voxel space Z^. 

Refen'ing now to Figs. 7 through 13, the apparatus of the present invention, will now be described. 
Notably, however, such apparatus while capable of carrying out the methods of the present invention, can 
nevertheless be used to psrfdrm other more general voxel-based and general computer data storage. 

18 
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accessing and processing operations. , . 

In Rg. 7. the voxel-based graphics system of the preferred embodiment, is graphically represented as 
comprising three processors, namely a Geometry Processor 4. a Frame Buffer Processor 3 and a Viewing 
Praeessor 5. a memory and addressing system, and a 2-D Frame Buffer 7. In general, the memory and 
addressing system 10 includes an addressing system 11. a Cubic Frame Buffer, and a voxel mulUple wnte 

'''"^R^a tte"»I«eH>a8ed graphics system of Fig. 7 is schematically illustrated, with the elements of the 
menury and addressing system 10 shown in greater detail. Specifically, the memory and addresang 
system 10 has ai overall modular stnjcture comprising a plurality of similar modules 13 indexed i(=o. i ... 
n-1 each central addressing unit 14. and several shared buses to be described hereinafter in greater detail. 

In the preferred embodiment, each k-th module 13 includes a k-th h)cal address unit 15. a k-th memory 
module 16. and voxel multiple write processing unit i7. Each k-th memory module is independer»«y 
accessible (i*. addressable) by the central addressing unit 14. as well as in cooperation with the 
con-esponding k-th local aWressing unit 15. 

in the preferred embodiment, each k-th memory module contains a plurality of internal memory cells 
each indexed (or addressed) by Internal memory cell indices (or addresses) i and j. where OfiSn-i and 
OSjin-1 Notably however, such a two^jimensional type memory unit is not necessary to carry out the 
method hereof, but couW be realized using a Bnear type memory array specified by a single index n 
ranging from OSnS(n-l)? for an n>-voxel Cubic Frame Buffer. 

In order to write In and read out voxel-based data with the plurality of n memory modules, a voxel data 
bus 18 is provided and can bo realized in several ways known to those skilled in the art. In additton. an 
intermodule voxel-data movement {i.e. shifting or skipping) mechanism 19 is also provided, being more fully 
illustrated in Rg. 12. Through each of the plurality of voxel-multiple write processing units 17. the operation 
and function of vrhfeh will be described in detail hereinafter. . . ...^ 

Referring now to Fig. 9. there is graphically represented in detail, the imernal elements of a single k-th 
module of the modular-stnictured memory and addressing unit iUustrated in Fig. 8. Also shown in Rg. 9, is 
(he central addressing unit 14 which, in the preferred embodiment hereof, provktes in a global manner, 
broadcasted addresses k". i'. and j' to each of the local addressing units of the n mpdujes^by way of a 
common bus for the purpose of transmitting (i.e. "broadcasting") those addresses k , i . j to the kwai 
30 addressing units, in a manner to be described more fuHy hereinaftor. 

Also shown in Rg. 9. are inteival and depth control parameters, e^.. dt and da. which are broadcasted 
to each of the local addressing units of the n modules, preferably over a common bus. In addition, voxel 
beam write^adMew control parameters are similarly "broadcasted" over a common bos. to each of the K- 

th tocal addressing units. , , ^ ^, 

3S In each k-th module, the tocal addressing unit IS includes several elements; namely: a local voxel depth 
measure (i.e. metric) computatton unit 20 for computing the voxel depth measure of the k-th voxel in a 
beam; a k-th module index fi.e. identity) storage unit 21; a k-th broadcasted address/module identity 
comparator unit 22 comparing the k' broadcast address with the k-th module index and generating a control 
enable: and comparalor unit 23 for comparing the -broadcasted" Interval and depth control parameters d, 
40 and d2. with the tocally^omputed voxel depth measure of the k-th voxel, and generating a control enable as 

well, to be provMed to the k-lh memory module. ^ , ... u 

Each k-th local addressing unit 15, is provkted with a hjgteal -OR" devtee 24, the inputs of wh.ch are 
the outputs of the comparator units 22 and 23. Depending on the mode of memory access of the memory 
and addressing system 10 (i.e. whether a single voxel or beam of voxels is to be read, wntten, or viewed), 
one or more than one hjgical "OR" device provides control enable signals to its corresponding memory 
module. For example, when only a single k-th voxel Is to be read or written, then only the k-th memory 
module will be enabled, tho k-th enabled memory modulo being the module where the broadcasted Indw 
k' equals the k-th module index (or identity). On the other hand, when a beam of voxels is to be read, 
written or viewed, then all of tho memory modules will be enabled, where the enabled memory modules are 
those which fall within the 1nten«l" specified by the -broadcasted" intenral control parameters d. and da. 
In either mode of memory access, the enabled memory module is enabled to read or wnte voxel(s) on the 
voxel data bus. 

Each local addressing unit IS also inchides a switching means for passing either 
fi) the broadcasted addresses l' and j'. ^ . , , u. 

ss m the addresses equated to the voxel depth measure d locally computed by local voxel depth 

'^''*2 comWnaton thereof as wHI be explained further hereinafter. In the preferred embedment shown 
in Hg. 9. the switching means is realized using a pair of multiplexers 25 and 26. the controls of which are 
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not shown. Notably however, the controls to these multiplexers 25 and 26 are specified, preferably, by 
"voxel beam direction** specifications selected by the observer/user on the basis of "spatial parameters** 
such as £x, £y. and tz, to be described in greater detail hereinafter. 

It is appropriate at this juncture, to point out that during voxel reading and writing modes of operation, 

s the memory module and local address units of each nrK)dule may participate in carrying out the specified 
operation, and that, in general, although subject to some exceptions, the voxel multiple write processing 
system is inactive. One exception Is during intermodule voxel-data movement, carried out by the inter- 
module voxel-data movement mechanism. There, tiie voxel-multiple write processing system is enabled to 
assist in performing the necessary operations of such an intermodule voxel-data permutation scheme, which 

fo will be descrit>ed in greater detail hereinafter. 

During the viewing mode of the memory and addressing system 10, ttie voxel multiple write processing 
system 12 is an active participant in such a projection-based process. 

Referring to Fig. 9, each k-th voxel multiple write processing unit 17 includes in the prefen^d 
embodiment a voxel value storage register 27. a voxel depth measure storage register 28. a local voxel 

75 depth competition unit 29. a transparency control unit 30. a switching unit 31 . and a logical "AND" device. 
Each k-th voxel value storage register 27. ts interconnected with each of its atfiacent neighboring voxel 
value storage registers by way of a voxehdata transfer patii. 

As illustrated in Ftg. 12. the intermodule connections amongst the voxel storage registers can ''skip" 
one or more voxel storage registers atong the voxel data transfer path. Such "skips** couM be realized fdr 

20 example, by hardwired connections. 

The output of each k-th voxel value storage register 27 is provided to two input ports of two elements, 
namely, the transparency control unit 30 and the switching unit 31. 

In addition, each k-ti^ transparency control unit 30 is provided with transparency control parameters 
transmitted (l.e., "broadcasted**) to all transparency control units, over a common bus 32. The nature of 

25 these transparency control parameters can foe of a wide variety, as can be the voxel value tiiemselves. For 
example, when using an 8-bit voxel representation, ttie voxel value can represent different colors. Thus, 
within this range of colors, any set of such colors can be defined as **transparent'*, or opaque (i.e.. non- 
transparent). It is the function of the transparency control parameter, tiien to define which cok>rs or 
combinations of colors are to be announced as "transparent" or *'opaque". rtotably. however, voxel values 

30 can represent physical or otiier properties of objects and scenes, other tiian color. For example, texture, 
material density and translucency ratios, of a small unit cube in the real scenes are Important. 

In order to carry out the projection method wherein tiie voxel (in a beam along in a conflict force 
storage retrieval direction) closest to an observer is determined, the method of the present invention 
involves tiie following procedure. For each voxel in a beam, the voxel depth measured along a specified 

35 storage/retrieval direction, is placed In its respective voxel depth measure storage register 28. Then each k- 
th transparency control unit determines whether the voxel value in its corresponding voxel value storage 
register is, by definition of the broadcasted transparency parameter, opaque, in addition, if tiie voxel deptii 
measure of tiie voxel in tiie corresponding voxel value storage register is determined to be within the range 
of deptii control (i.e.. in Uie "undipped** range) specified by d. and da. tiien a deptii control enable is 

40 provided to the input of the togical "AND" device. If tiie voxel value In tiie k4h voxel value storage register 
is botii "opaque" and within tiie deptii interval specified by d, and 62 control parameters, then tiie and 
device generates a deptii competition enable signal to tiie k-tti corresponding voxel deptii competition unit, 
ttieretiy allowing it to conduct a local **deptii value competition process** to be described below. Notably, 
hovyever. tiie "AND** device can be substituted witti an "OR" device or eliminated altogetiier ttius allowing 

45 k)cal deptii competition processes In each of the n modules, regardless of ti'ansparency and deptii interval 
parameters. This might be useful on reading/writing operations. Notably, this ••winning** voxel value 
represent tiie "opaque" voxel value which closet to tiie observer viewing along the voxel beam parallel to 
tiie user*specified storage/retrieval direction. In addition, tiiis process provides great flexibility in detining 
whtoh voxel values are opaque, and tiius are allowed to compete locally prior to passage into tiie 2-D frame 

50 buffer, where tiie voxel value can be processed "on tiie fly" wittiln a suitable shading unit. 

Returning now to Ftg. 11, tiie internal elements of each voxel multiple write processing unit 17 of ttie 
preferred embodiment, is illusti^ated in greater detail. 

As illustrated, the k-tti local voxel depth computation unit 20 provkies to the k-thvoxel deptii measure 
storage register 28, ttie computed deptii measure dk. which as shown. Is represented by a plurality of (e.g. 

55 nine) binary bits, bg, b?-. bi...bo, for ttie case of n » 512. The depth measure storage register 28 comprises 
a plurality of storage cells, each cell containing a bi bit which is provided to a different stage of the k-th 
voxel deptii competition unit 29. Each t-tti stage in the voxel deptii competition unit comprises conventional 
togk: circuitry (e.g., logic elements A. B. and C) which in effect compares ttie bi btt witti ttie corresponding 
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b) bits of all other voxel depth measures provided to the k-th depth competition unH by way of the commoti 
depth bus. as described hereinbefore. 

Tha operation of the voxel depth competition unit 29 is as follows with the (i=8)th stage of each 
-enabled" voxel depth competition unit being active while all other stages are inactive, the most significant 

5 binary bit digit of each voxel depth measure is compared using the logic circuitry of that stage, and the 
group of voxels having the predetermined extreme value (e.g.. largest value) are determined. For each b. bit 
(in each depth competition unit 29) which was determined to have the predetermined extreme value, me 
next stage of these depth competition units are enabled for a comparison of the next most significant bit 
value the above comparison, determining and group formation steps are repeated for successive stages of 

to the voxel depth competition unit 29. on the basis of consecutive significant binary digits, until only one 
surviving voxel depth measure remains in the group. The voxel depth competition unit 29 corresponding to 
the surviving voxel, generates a voxel pass control provided to the switch unit 31 to allow that voxel in fte 
corresponding voxel value storage register, to be passed into 2-0 frame buffer. (FB) 7 as descnbed 
hereinbefore. The surviving voxel depth measure corresponds to the voxel which is closest to the viewnng 

IS along the specified storage/retrieval direction. ^ ^ 

Returning to Fig. 13 in particular, an addition schematic representation is shown. lOusiratfng only those 
portions of central addressing unit local address units, and memory modules which funcBoo in ttie memory 
storage mapping and demapping operations of the present invention. As discussed hereinbefore, pairs of 
coordinate value sums and differences are required to compute the local addresses using me demapping 

20 functions as discussed hereinbefore. However, since such coordinate value sum and differences are 
required by each local voxel depth computation unit, they are in me prefened embodiment, comjwted by 
the central addressing unit 1 4 and "broadcasted" to all local addressing units over a common bus 32. 

As illustrated in Fig. 13.. me demapping functions (i.e.. voxel depm measure fonnulas) can be realaed 
using multipliers, adders, subtractors and switches. For me case of six viewing directions, a depm measure 

25 function is realized in Rg. 13. partly by me central addressing unit U and partly by me k-me tocal 
addressing unit 15. However, modifications of such realizations are expected to occur. 

In addition while not shown in Fig. 13. the binary skewing function S (x. y, 2) for me case of six conflict- 
free storage/retrieval directions, could typically be realized in the central processing unit 14 using again, 
devices to carry out multiplication, addition, subtraction and other arimmetlc functions, and thereafter be 

30 "broadcasted" to all in local addresses, over a common addressing bus 32. 

While me particular embodiments shown and discussed hereinabove have proven to be useful m maiiy 
applcations. further modifications of the present invention hereindisclosad will occur to persons skilled m 
me art to which me present invention pertains, and all such modifications are deemed to be within me 
scope and spirit of me present invention defined by the appended claims. 



Clatms 



1 A method of storing and retrieving one or more voxels of a beam disposed along one or more of a 
40 plurality of storage/retrieval directions of a 3-0 matrix array in 3-0 discrete voxel space, said 3-0 discrete 
voxel ^)ace being specified in terms of x. y and z coordinate directkm. said method comprising the steps 
of' 

(a) storing said voxels into memory storage space by mapping said voxels along one of said plurality 
of storage^retrieval directions into a plurality of n independanOy accessible memory modules m said 

43 memory storage space, each said memory module being indexed by a memory module index and having 
internal memory cell addresses, said n being an integer, said mapping being carried out in accordance wrth 
a skewing function expressed in terms of said x. y and z coordinate directions and saU Integer n; md 

(b) retrieving said voxels along one of said plurality of storage/retrieval direcfions in said 3^) voxel 
space, by demapping said one or more voxels from said memory storage space into saM 3-D voxel space 

90 using spatial parameters and said integer n. ^ jm. r „ 

2. The memod of claim i. wherein step (a), said mapping is carried out in accordance wrth a linear 
skewing function expressed in terms of said x, y and z coordinate directions and said integer n. 

3. A memod of claim 2. wherein step (a) comprises: 

0) mapping sakl voxels along said beam parallel to sakl storage/retrieval directton. into a plurality of n 
56 independently accessible memory modules being indexed vrim memory module indices k=0. 1. ••• "-J. 
each said k-th memory modute conteining memory cells, each said memory cell being specWied by 
intemat module indices, i and j. said mapping being carried out in accordance wHh a linear skewing function 
of the general form 
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k ■ (ax + by + cz) mod n 

where x. y and 2 are the coordinate values of each said voxel in said 3-D discrete voxel space, where i =x 
and jay. and where a, b and c are integer coefficients and n is a prime relative to said integer coefficients 
a. b and c. and linear combinations of said integer coeffidonls a, b and c the number of storage/retrieval 
5 directions in said n^ matrix anray in said 3-D discrete voxel space being based upon the values of said 
integers a, b. c and n; and wherein step (b) comprises: 

(i) for each said voxel along said beann parallel to said storage/retrieval direction, determining a voxel depth 
measure representative of the depth of each ssdd voxel in said n^ matrix array in said 3-D discrete voxel 
space, said determination of said voxel depth measure being determined for each said voxel on the basis of 

10 one or more of said x, y and z coordinate values of said voxel, said integer coefficients a. b and c. and said 
module index k of said memory module into which said voxel has been mapped in step (a), 

(ii) for each voxel along said beam parallel to said storage/retrieval direction, detennining said memory cell 
indices I and i on the basis of said voxel depth measure and said x. y and 2 coordinate values 
corresponding to each said voxel, and 

75 C«0 retrieving said voxels along a specified storage/retrieval direction, using said memory ceil incfices i and i 
determined in step <ii) and said module index k determined in step (a) 0)- 

4. The method according to claim 3. wherein each said integer coefficients a. b and c equals 1. so as to 
provide six storage/retrieval directions corresponding to the principal axes ±x. ±y. tz of said n^ matrix array. 

5. The method according to claim 4. wherein step (b) (i) comprises detennining for each sakl voxel. 
20 said voxel depth measures d*x and d.x for said voxels along sakJ beams parallel to the x principal axis, said 

voxel depth measure d*y and d^, for said voxels along said beams parallel to the y principal axis, and said 
voxel depth measure d*, and d., for said voxels along said beams parallel to said z principal axis, being 
determined, respectively, according to the functions 
d*x = (k-(y ♦ z))modn 
25 d^ = n • d*K- 1 

d*y = (k • (x + 2)) mod n 
d-y =* n - d4>y • 1 

= (k - (X ♦ y)) mod n 
d^ « n - d*j- 1 

30 a The method according to claim 3. wherein said integer coefficient a equals 5, said integer coefficient 

b equals 2. and said integer coefficient c equals 1. so as to provide twenty-six storage/retrieval directions of 

sas6 n^ matrix array in said 3-D discrete voxel space. 

7. A method of determining which voxel in a plurality of voxels along a beam parallel to a 

storage/Mrieva! direction of a matrbc array in 3-D discrete voxel space, is closest to an observer viewing 
35 atong s^d beam, saki 3-D discrete voxel space being specified In terms of x. y and 2 coordinate directtons» 

said method comprising the sequence of steps of. 

(a) mapping said voxels atong said beam parallel to said storage/retrieval direction Into a plurality of n 

independently accessible memory modules being indexed with module indices k=0, 1 n-i. each said k- 

th memory module containing memory cells, each said memory cell being specified by internal module 

40 indices i and i. sakJ mapping being carried out in accordance with a linear skewing function of the general 
form 

k 3 (ax by +cz) mod n 

where x, y and 2 are the coordinate values of each said voxel in said 3-D discrete voxel space, where I =x 
and j"y. and where a. b and c are integer coefficients, the number of storage/retrieval directions in said 
45 matrix array in said 3-D discrete voxel space being based upon the values of said Integers a. b. c and n; 

(b) for each said voxel in said string of voxels along sakj beam parallel to saW storage/retrieval 
direction, determining a voxel depth measure representative of the depth of each sakl voxel in said matrix 
array In said 3-0 discrete voxel space, said determination of said voxel depth measure being determined for 
each said voxel on the basis of one or more of said x. y and 2 coordinate values of said voxel, said integer 

50 coefficients a, b. and c. and sakJ module index of sakI memory module into whtoh sakj voxel has been 
mapped in step (a); and 

(c) comparing said voxel depth measures determined in step (b). so as to determine which voxel 
depth measure a predetermined extreme value, and thereby representative of the con'esponding voxel 
being closest to said obsen/er viewing along said beam parallel to said storage/retrieval direction. 

S5 8. The method of claim 7, wherein said voxel is further specHied by a transparency parameter which is 
either a transpaient or a non-transparent voxel value, and wherein said method further comprises prior to 
step (c). 

determining said voxels having a non^transparent voxel value, and wherein step (c) compnses 
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comparing only sard voxel depth .measures determined in step (b) which con-espond to said voxels 
determined above as having a non-transparent value, so as to determine which voxel having a non- 
transparent voxel value, has a predetermined extreme voxel depth measure value, and thereby representa- 
tive of said voxel being closest to said observer viewing along said beam parallel to said sft>rage/retrievat 
s direction. 

9. The method of claim 10. wherein said method comprises, prior to step (c). 
determining said non-transparent voxels along said beam, which have corresponding voxel depth measures 
within at least one predetermined depth measure range 
wherein step (c) comprises 

10 comparing only said voxel depth measures detemiined in step (b) which con-espond to said non-transparent 
voxels detemilned above as having corresponding voxel depth measures within said predetermined depth 
measure range, so as to determine which voxel within said predetermined depth measure range, has a 
predetermined extreme voxel depth measure, and thereby representative of said voxel being closest to said 
observer viewing along said beanft parallel to said storage/retrieval direction. 

IS 10. The method of claim 7, wherein each said voxel depth measure is represented by a plurality of 
binary digits, and wherein said comparison of voxel depth measures comprises 

(i) comparing the most significant binary digit of each said voxel depth measure. 

(ii) determining which voxels have the predetermined extreme voxel depth measure, on the basis of 
comparing said most significant binary digits, 

20 (iii) forming a group of voxels from the determination performed in step (ii). 

(iv) comparing the next consecutive most significant binary digit's of voxel depth measures of said 
voxels in said voxel group formed in step (iii). 

(v) determining which voxels in said voxel group have the predetemnined extreme voxel depth 
measure, on the basis of comparing said next consecutive largest significant binary digits, and 

26 (vi) repeatedly perfomrting said comparing, determining and group formation steps on the basis of 

consecutive significant binary digits, until only one sun/ivmg voxel depth measure remains in said group, 
said sun/lving voxel depth measure conresponding to said voxel being the closest to said obsenrer viewing 
along said storage/retrieval direction. 

11. The method of claim 10, which further comprises. 

30 for a plurality of beams parallel to said storage/retrieval direction, each said beam having a plurality of 
voxels therealong. 

(1) mapfwng said plurality of voxels along each said beam, into said plurality of n independently 
accessible memory modules according to said linear skewing function. 

(2) determining a voxel depth measure for each said non-transparent voxel along each said beam. 

3S (3) for each said voxel along each said beam, comparing said voxel depth measures corresponding 

to said non-transparent voxels, so as to determine which voxel depth measure is a predetermined extreme 
value, and thereby representative of the corresponding non-transparent voxel being closest to said observer 
viewing along said beam parallel to said storage/retrieval direction, and 

(4) fonning a two-dimensional image from said plurality of non-transparent voxel values detennined in 

40 step (3) above. 

12. A method of mapping the voxels of an n^ matrix array in 3-0 discrete voxel space, into a plurality of 
n independently accessible memory modules, said 3-D discrete voxel space being specified in temns of x, y 
and 2 coordinate directions. 

said method comprising: 

45 (a) defining a plurality of 3-D voxel subspaces sequentially numbered module n. for K = 0, 1, ... n-1. 

and grouping together each said K-th 3-D voxel subspace. each said K-th 3-D voxel subspace containing a 
set of voxels having x. y and z coordinate values; 

(b) providing a plurality of independently accessible nxn memory modules indexed K=0. 1. ... n-1, 
each said K-th memory module containing r? memory cells, each memory cell being specified by inlemal 

50 module indices, i. j; 

(c) assigning each grouping of saW K-th 3-D voxel subspaces to said K-th memory module; and 

(d) for each said voxel having said x, y and z coordinates, mapping from said 3-D voxel space, into 
said K-th memory module according to 

K = (5x + 2y ■> z) mod n 
85 and Into the 0. j)-th memory cell according to x^i and y sj. 

13. Apparatus for simultaneously storing and retrieving a beam of voxels disposed parallel to a 
storage/retrieval direction of an n^ matrix an-ay in 3-D discrete voxel space, said 3-D discrete voxel space 
being spedfidd in terms of x, y and z coordinate directions, said apparatus comprising: 
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storing means for storing said voxels into memory storage space, by mapping said voxels along one of said 
plurality of storage/retrieval directions, into a plurality of n independently accessible memory modules in 
said memory storage space, each said memory module being indexed by a memory module index and 
having internal memory cell addresses, said n being an integer, said mapping being carried out in 
5 accordance with a skewirtg function expressed in terms of said x. y and z coordinate directions and said 
integer n; and 

retrieving means for retrieving said voxels along one of said plurality of storage/retrieval directions in said 3- 
D voxel space, by demapping said one or more voxels from said memory storage space into said 3*D voxel 
space using spatial parameters and said integer n. 
10 14. Apparatus of claim 13. wherein said storage and said retrieving means comprises: 

(a) a plurality of n modules being indexed with module indices k = 0. 1. ... n-1, each said Mh module 
including a k-th memory module, a k*th Ickcal address unit for independently accessing said k-th memory 
module, each said memory ceil of each said k-th memory module being specified by internal module 
indices i and j, each said k-th local address unit including a \ocb\ voxel depth computation unit fbr 

15 computing a voxel depth measure of each voxel residing in said k-th memory module, said voxel depth 
measure of each voxel being representative of the depth of said voxel in said n^ matrix array, and being 
determined on the basis of one or more of said x, y and z coordinate values of satd voxel, said integer 
coefficients a, b and c, and said module index k of said memory module into which said voxel is mapped, 
said local addressing unit further including means for determining for each said voxel, said memory cell 

20 indices i and j on the basis of said voxel depth measure and said x. y and z coordinate values 
corresponding to each said voxel, and memory cell index selection means for selecting memory cell indices 
I and j from global memory cell indices i' and j " and said voxel depth measures ak)ng said storage/retrieval 
directions. 

(b) a central control unit Including a central addressing unit for generating module indices k'»0, 1, ... 
25 n-1. and gtobal memory cell indices i" and j", said central addressing unit further Including means fbr 

canrying out the mapping of said voxels along said beams parallel to said storage/retrieval direction, into 
said plurality of n independently accessible memory modules, said mapping being carried out in accor- 
dance with a linear skewing function of the general form 
k a (ax + by + cz) mod n 

30 where x, y and z are the coordinate values of each said voxel in said 3-0 discrete voxel space, where i^x 
and jsy. and where a, b and c are integer coefficients and n is a prime relative to said integer coefficients 
a, b and c, or to linear combinations of said Integer coefficients a. b. and c of the number of 
storage/retrieval directions in said n^ matrix array being based upon the values of sakJ integers a. b. c and 
n. 

as (c) means for transferring global memory cell indices i" and j" from said central addres^ng unit to 

sakj plurality of n local address units, and 

(d) means for transferring said beam of voxels, into or from said n memory modules. 
15. The apparatus of claim 14, wherein each said k-th module further comprises 
a voxel multiple write processing unit for processing voxels accessed from said k-th memory module, 
40 a k-th voxel value storage register for storing the k-th voxel value of the k-th voxel atong said beam parallel 
to said storage/retrieval direction. 

a k-th voxel depth measure register for storing said voxel depth measure of said k-th voxel along said beam 
parallel to said storage/retrieval direction, 

a k-th voxel depth competition unit for locally comparing the value of said voxel depth measure of said k-th 
45 voxel with all other voxels of said beam parallel to said storage/retrieval direction, and determining whtoh 

said voxel depth measure has a predetermined extreme value and generating a voxel pass control signal 

corresponding to said voxel depth measure having said predetermined extreme value, and 

writing means for writing into a two*dimensional frame buffer, said voxel value stored in said k«th voxel 

value storage register corresponding to said voxel pass control signal. 
50 16. Apparatus of claim 15. which further includes means for providing transparency control parameters 

to said voxel multiple write processing units, each said k-th voxel multiple write processing unit further 

including 

a k-th voxel transparency control unit for receiving said transparency control parameters and for controlling 
the local comparison of said voxel depth measures within said k-th voxel depth competition unit saki voxel 
55 depth competition control being earned out on the basis of said voxel value in said k-th voxel vahje register 
and said transparency control parameters. 

17. The apparatus of claim 14. wherein each said k-th local address unit includes a tocat comparison 
means for comparing the voxel depth measure of each said voxel with at least one preselected voxel depth 
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measure range, and controlling the local comparison of said voxel depth measures within said k-th voxel 
depth competition unit, said voxel depth competition control being carried out on the Iwsis of said voxel 
tjeing within said preselected voxel depth measure range. 

18. Apparatus of claim 15. wherein each said voxel depth measure is represented by a plurality of 
5 binary bits, and wherein said local voxel depth competition unit comprises means for 

(i) comparing the most significant binary bits of each said voxel depth measure. 
Oi) determining which voxels have a predetermined extreme voxel depth measure, on the basis of 
comparing said most significant binary bits. 

(lit) forming a group of voxels from the determination performed in step (ii). 
10 (iv) comparing the next consecutive most significant binary bits of voxel depth measures of said 

voxels in said voxel group formed in step (iii). 

(V) determining which voxels in said voxel group have the predetermined extreme voxel deptn 
measure, on the basis of comparing said next consecutive most significant binary bits. 

19. The apparatus of claim 15. which further includes intermodule voxel value moving means for moving 
15 voxel values stored in said k-th voxel value storage registers, to other voxel value storage registers. 

20. A 3-D voxel-based graphics system comprising: 

a memory and addressing system for storing and retrieving a 3-0 voxel-based image along a plurality of 
storage/retrieval directions in 3-0 discrete voxel space, said memory and addressing system Including a 
plurality of modules each said module having a memory module, and a local address unit for independently 
20 addressing said memory module; ^ * , 

a central addressing unit for transmitting to each said module, broadcasted addresses and control 

parameters; • a 

means for transferring said broadcasted addresses and control parameters to each said module; ana 
means for transferring said voxel-based image into or from said n memory modules. 
25 21. The 3-0 voxel-based graphics system of claim 20. wherein said system further comprises 
at least one voxeMased processor for accessing said n memory modules. 

22. The 3-D voxel-based graphics system of claim 21, wherein said voxel-based processor is selected 
from the group consisting of a 3-D geometry processor, a 3-D frame buffer processor, and 3-D viewing 
processor. 
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© The present invention concerns a method of and 
apparatus for storing to and accessing from a mem- 
ory device, one or more voxels of a t^eam disposed 
along one or more of a plurality of storage/retrieval 
directions of a 3-D matrix an^ay and 3-D discrete 
voxel space. Voxels are stored into memory space 
by mapping the voxels along one of the plurality of 
storage/retrieval directions, into a plurality of inde- 
pendently successful memory modules in memory 
storage space. Each memory module is indexed by 
a memory module index and has internal memory 
cell addresses. The mapping is carried out in accor- 
dance with a linear skewing function which ex- 
pressed in terms of x. y, and z coordinate directions 
and integer n. The voxels can be retrieved along one 
or more of a plurality of storage/retrieval directions 
by de-mapping one or more of the voxels from the 
memory storage space into 3-D voxel space using 
spacial parameters and integer n. The demapping 
operations of the present method are also used in 
order to determine the voxel depth measures of 
voxels alCHiQ one of the storage/retrieval directions. 
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